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Photon correlations and cross-correlations of light scattered by a regular structure of strongly 
driven atoms are investigated. At strong driving, the scattered light separates into distinct spectral 
bands, such that each band can be treated as independent, thus extending the set of observables. 
We focus on second-order intensity-intensity correlation functions in two- and multi-atom systems. 
We demonstrate that for a single two-photon detector as, e.g., in lithography, increasing the driving 
field intensity leads to an increased spatial resolution of the second-order two-atom interference 
pattern. We show that the cross-correlations between photons emitted in the spectral sidebands 
violate Cauchy-Schwartz inequalities, and that their emission ordering cannot be predicted. Finally, 
the results are generalized for multi-particle structures, where we find results different from those 
in a Dicke-type sample. 

PACS numbers: 42.50.Hz, 42.50.St, 42.50.Ct 



I. INTRODUCTION 



Throughout the history of quantum physics, first-order correlations have played a major role in discussing the 
foundations of the underlying theory. One of the most famous model system is Young's double-slit, which despite its 
simplicity allows to explore fundamental questions such as complementarity and uncertainty relations [ll [3, la, 3l ■ A 
modern realization of Young's experiment involves two atoms scattering near-resonant laser light [E,@,0,HHI3|- The 
light is thus scattered by the simplest form of a regular structure, which gives rise to interference phenomena in the 
scattered light, because different indistinguishable pathways connect source and detector Analogous interference 
is also possible with single particles, demonstrating pathway interference in the energy-time domain (T^ |. 

Starting with the experiments by Hanbury-Brown and Twiss [l3|] involving the measurement on a thermal field 
and the corresponding measurement with anti-bunched resonance fluorescence light [T3 | , the second-order correlation 
function gained considerable attention. This correlation function measures intensity correlations and thus provides 
information on the fluctuation of the fields [IB]] . 

In general, spatial interference in the different correlation functions can only be expected under certain experimental 
conditions. For example, it is possible to have interference in the second-order correlation function under conditions 
of no interference of the first-order correlation function Q ■ In particular the first order interference in light scattered 
of re gula r structures of atoms is typically restricted to low incident light intensity and vanishes at strong driving 0, 
H, l& llCt [H, H?], UK]. For example, in a two-particle system in the strong field limit, the collective dressed states 
are uniformly populated, and thus the interference fringe visibility is zero. Since applications often rely on these 
interference effects [H, [It], HH Ojj] , this restricts the potential implementation, where properties such as coherence, 
increased resolution, high signal-to-noise ratios or a rapid coherent system evolution are among the desirable properties 
of the scattered light. 

More fundamentally, different lines of the frequency spectrum of the scattered light become separated in the strong- 
field limit, just as in the Mollow resonance fluorescence spectrum of a single strongly-driven two-level atom. We have 
shown recently that this separation can be employed to recover full first-order interference fringe visibility in the 
strong- field limit, and to gain a clearer interpretation and an extended set of observables to describe the scattered 
light [18l | . The basic idea was to facilitate a frequency-dependent electromagnetic bath, which can be realized, for 
example, using cavities or photonic crystals [20I. |2lj. The focus of this work, however, was on first-order correlation 
functions. 

Here, we discuss spatial second-order quantum interference effects in strong driving fields. As in the strong- field 
limit the spectral lines are well-separated, we define observables for each of the spectral lines separately. This allows 
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for additional observables such as cross-correlations between different the spectral lines. First, we investigate the 
spatial dependence of the second-order correlation functions for a strongly driven atomic pair in free space and focus 
on the case of a single two-photon detector registering photons at the driving laser frequency u>l, as, for example, 
in lithography with a medium sensitive to two-photon exposure. We show that in this setup, the spatial resolution 
of the central strong-field second-order interference pattern can be increased by a factor of two as compared to the 
corresponding weak-field pattern simply by increasing the driving field intensity. This allows to create structures with 
high spatial resolution and signal intensity. Next, quantum cross-correlations between photons emitted in different 
spectral bands are investigated. In particular, we show that the spatial Cauchy-Schwarz inequalities (CSI) are violated 
for photons emitted into the sideband spectral lines over a wide range of detector positions. Also, it is impossible to 
predict the temporal ordering of two photons emitted in different sidebands. Our scheme can be realized in a wide 
range of systems, and can also be employed to analyze the structure of the scatterers. Finally, we generalize our 
results to the case of a linear chain of N atoms, where the results are found to be different from results in Dicke-type 
samples. 

II. ANALYTICAL TREATMENT 

We first concentrate on a driven pair of distinguishable non-overlapping two-state emitters, labeled a and 6, and 
generalize the discussion to a linear structure of N > 2 atoms in section|Vl The two atoms both have atomic transition 
frequency loq, are located at positions r a , fb, and are separated by r a \,. The electronic bare ground state of atom 
j € {a, 6} is labeled by the corresponding excited state is \2}j. The external driving field with frequency ll>l = cHl 
and wave vector k^ is aligned such that k^ ■ f a b = (see Fig. [T]). The particles are spontaneously damped via the 
interaction with the surrounding electromagnetic field (EMF) reservoir. Our aim is to investigate the spatial photon 
correlations of the light scattered by the atoms in the intense driving field limit. 

Assuming the electric dipole and rotating wave approximations, the laser-dressed atomic system is described by the 
Hamiltonian H = Hq + Hj, with 

H Q = ^h(uj k - LLi L )a[a k + ^ h&Rzj, (1) 

k je{a,b} 

H T =iJ2 E • & \ a l ( R *j^r - R 2i sin2 & + Ri2 cos 2 o) e -«) ^ - H.c.) . (2) 

Here, Hq is the Hamiltonian of the free EMF and the free dressed atomic subsystems, while Hi describes the interaction 
of the laser-dressed atoms with the EMF. a k and a k are the radiation field annihilation and creation operators obeying 
the commutation relations [ak, n!/] = 8 kk > , and [ak, a k'}— [ a L a l'] = 0- The atomic operators — \a)jj(l3\ describe 
the transitions between the dressed states \0)j and \a)j in atom j for a ^ (3 and dressed-state populations for a = (3, 




FIG. 1: A linear structure of two-state emitters, where neighboring atoms are separated by the distance vectors r a i. The 
whole structure is driven by a resonant strong external laser field with wave vector &l. The two detectors are used to measure 
correlations among the photons emitted in observation directions _Ri,i?2, which are also characterized by the angles ai,«2 
between r a t, and R\,R2- The line of atoms and the two detectors are located in a plane, and the driving field propagates 
perpendicular to this plane. 
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and satisfy the commutation relation 



l^lR ( l>] = Sji S aa 'R^,-S 0la R% . (3) 



V» 1 V/3'J ~~ 



>C0 _ x . R 0) 



a/3' " V/3 



(7) 

The operators -R^g can be represented through the bare state operators via the transformations 

= sindfyj +cos0|l) 3 - , \2),j = cos0|2)j - sin^l^ , (4) 

where the mixing angle 9 is given by cot 29 — A/ (217). The laser field detuning is A = ujq — lol, and the Rabi frequency 
is defined by 2f2 = (d ■ E£)Jh. Here, El is the electric laser field strength, and d = d a — db is the transition dipole 
matrix element. We further define the population inversion operators R Z j = \2)jj(2\ — and the generalized 

Rabi frequency f2 = [Q 2 + (A/2) 2 ] 1 / 2 . The dressed state transition frequencies are ujl, w± =Wi±2f2. The two-particle 
spontaneous decay and the vacuum-mediated collective interactions are given by the frequency-dependent expression 

Ijl M = 7M [Xjl M + iflji (w)] ■ (5) 

Independent of the atom- vacuum coupling, the collective parameters ^ and fijj (j I) tend to zero in the large- 
distance case Tji — + 00 which corresponds to the absence of coupling among the emitters. In the small-distance case 
Tji — > 0, the parameter Xjl tends to unity, while "fflji tends to the static dipole-dipole interaction potential. 

In the intense field limit Cl 3> N-f, different lines of the spectrum are well-separated. This allows us to investigate 
the spatial dependence of the second - order correlation functions for each of the spectral band separately. In this 
case, it follows from the interaction Hamiltonian Eq. @ that the operators 

R zj sin29/2, cos 2 9 and R$ sin 2 9 (6) 

can be considered as the sources of the mth spectral field components {to £ C, R, L}. Here, C stands for the central 
spectral band emitted at frequency lul, and R and L are for the right and left sidebands emitted at frequencies Wi + 2f2 
and ojl — 2Q, respectively [22|. In the following, we use this decomposition to investigate the spatial properties of the 
photon statistics and cross-correlations between photons emitted from different spectral lines. 



III. TWO-PARTICLE QUANTUM DYNAMICS 

In order to be able to derive analytic expressions for the required expectation values, we first investigate the two- 
atom quantum dynamics in a strongly driven laser field. We introduce, for this reason, the two-atom collective dressed 
states as 

|* e ) = |2 a ,2 6 ), |*, (0) ) = -^(|2 a> U>±|2 6> l B )), and|* fl ) = |i a> U). (7) 

The variables a a p — (|^ r ct )(^ r | a|) are expectation values of the corresponding transition (a ^ (3) and population 
operators (a — (3) ({a, f3} £ {e,s,a,g}). Using Eqs. (|1I2[) . one can obtain the equations of motion for the dressed- 
state variable of interest, i.e., x — 2(a ee — cr gg ), y — a ss — <j aa , and z — a ee + a gg — a ss — a aa , as (23[: 

x(t) = ~2^x + AC i ab ) y + ^-\ (8a) 

i(t) = -&x 2(c$ + C (+) )y + (8b) 

z(t) = 2^x + 4^+ ) y-4^+h. (8c) 
The coefficients in these equations are defined as 

£ (±) = 7 (w_) sin 4 9 ± 7(0;+) cos 4 9 , (9a) 

Qb = l(u-)Xab(u-) sin 4 9 ± 7(w + )xah(w+) cos 4 9 , (9b) 

c ab = Tk)[l " XabM] sm 2 29 . (9c) 

Simple analytical expressions for the two-atom steady-state quantum dynamics can be obtained if the two-particle 
interaction is mediated by the usual vacuum modes of the environmental electromagnetic reservoir, i.e. 7 (w_) = 
7 (w+) = 7 (^i) = 7. Then, one obtains that 

x = 2£i-)/{M i y = Q, and z= K (_) /^ (+) ] 2 - (10) 

In this case, the diagonal atomic dynamics is independent of the inter-atomic separation providing that Q 3> 
{ 7 r2 a b, -/V7}. In particular for 9 = 7r/4, i.e., exact resonance, we find that x = y = z — 0, which means that 
the single-atom dressed states are equally populated. 



4 




FIG. 2: (a) Second-order side-band photon correlations g^l^^iRi, R2) and (b) cross-correlations 9^ R ^ RL ){Ri, R2) as function 
of the detector positions ai, 02- The interatomic distance is r a b/A = 5, and the number of atoms is TV = 2. 



IV. SECOND-ORDER CORRELATION FUNCTIONS AND TWO-PARTICLE SIDE-BAND PHOTON 

CORRELATIONS 



We now turn to the second-order correlation function of the steady-state resonance fluorescence emitted in the 
three spectral bands. The coherence properties of an electromagnetic field, at space-point R, can be evaluated with 
the help of the second-order coherence functions: 

(2)i Z B\ (gmfe jjjj °" ( f + T ' ^ + T ' am ^' jji) 
0,1 ' (T, iti, il 2 j = ; =; =; ; 5 =; , (11) 

(a^(«, J Ri)a m («, J R 1 ))(a+(i )J R 2 )a n (i,i? 2 )) 

where a+ (a„) are the photon creation (annihilation) operator for modes n S {C, L, i?}. The quantity g m n{ T ) can be 
interpreted as a measure for the probability for detecting one photon emitted in mode m and another photon emitted 
in mode n with delay r. From now on, all correlation functions are evaluated for r = 0, and, for notational simplicity, 
we drop the variable r. We further define two Cauchy-Schwarz parameters 



v ( f> 6 \ _ 9 L l( r ^ R 2)9 RR (Ri, R2) ( s g v _ g LL (Ri,R2)g RR (Ri,R2) ri9 , 

IrPLCR, P„ 2 [„(_ 2 CB, P„ 2 



which relate the correlation between photons emitted into individual modes to the cross-correlation between photons 
emitted into two different modes. If \l < 1 or \R < 1> the respective Cauchy-Schwarz inequalities are violated [24j . 
We now turn to detection in the far-zone limit, and specialize to a resonant driving field (9 = 7r/4). Then the field 

operators (a n ) entering in £ and corresponding, respectively, to the three spectral modes of the Mollow spectrum 
can be represented via the atomic operators given in Eq. ©. With the help of Eqs. (fTl fTU)) . and in Born-Markov and 
secular approximation, one arrives at the following expressions for the two-particle photon-correlations: 

gQ C (R\,R-2) — 1 + cos Si cos 82 , (13a) 

9 { lI(Ri,R2) = g$i(Ri,R*) = ^[1 + 008(5! - 82)) , (13b) 

gf R {Ri,R2) = g^l(Ri,R 2 ) = \^ + cosft + S 2 )} , (13c) 

gPxi&iM) = 9xc(RuR2) = 1 for X e {L, R} . (13d) 

Here we used that in the resonant strong-field limit one has (R Z j) s — (R% : )s = (^21 )« — and (R Z iR Z j) s = 

(-R^i-Rif )s = 0, where i ^ j G {a, 6} denote one of the two atoms. The parameters 8d (d 6 {1, 2}) are related to the 
detection angles ad by Sd = kr a {> cos ay, where we have assumed a planar geometry of atomic sample and observation 
detections. 

Already the structure of Eqs. (|13p reveals interesting insight in the nature of the light emitted in the various 

(2) ~* ~* (2) — ' -* 

spectral bands. From the properties of the cos function it follows that < g cc [Ri, i? 2 ) < 2, and < g L [(Ri, i? 2 ) = 
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9rr(Rii R2) — 1) an( i 1 — 9lr(Ri' Rz) — 5fli(^i'^2) < 2. The cross-correlations involving the central spectral 
band are always unity and do not exhibit a dependence on the detection position. This immediately fixes the possible 
photon statistics of the emitted light. We observe a tendency to have pairs of photons emitted in the central band 
or pairs where one of the photons is in cither sideband (left and right). In general, for these correlation functions, 
values below unity denote sub-poissonian light statistics, a value of units indicates poissonian statistics, and values 
above unity stand for super-poissonian light statistics. 

(2) (2) 

An example for the second-order sideband photon correlation functions g LL and g RR as a function of the detection 

angles <5i, 82 is shown in Fig. 03a). The spatial interference fringes of the second-order correlation function are clearly 

(2) (2) 

visible, and the range of possible values for g LL and g RR shows that both poissonian and sub-poissonian photon 

statistics can be generated at the sideband frequencies. Fig.^b) shows the corresponding results for the second-order 

(2) (2) 

cross-correlations g LR and g RL - Here, the oscillatory structure is different, and values from 1 to 2 are obtained for 
the cross-correlations. Therefore poissonian as well as super-poissonian light statistics can be generated in the cross- 
correlations. In general, sub-Poissonian and Poissonian light statistics can be generated in all three spectral lines, and 

(2) 

super-Poissonian light statistics is possible in the central spectral band second-order correlation function g c ^, and in 
the cross-correlations, see, e.g., Eq. (ITcH) . 

We now discuss the second-order spatial interference of the light emitted in the central frequency band around 
ljl for a specific detection scheme. We assume identical detection angles 8\ — 82 = 8 and detection positions 
R\ = i?2 = R, which corresponds, e.g., to a medium sensitive to two-photon exposure fl9j |. In the following, we 
discuss the interference fringe resolution of this setup, which is of relevance to applications in lithography, where 
the general aim is to create structures as small as possible. In the strong- field limit fi/7 ^> 1, from Eq. (fT5|) we 

find g^ c {R) = 1 + cos 2 8. In the weak field case fi/7 < 1 without spectral band separation, however, one finds 
g^^R) — [s/(s + cos 8)] 2 with s = l + 2(fi/7) 2 [10| . Comparing these two results, one finds that simply increasing the 
driving field strength effectively doubles the spatial fringe resolution in this setup around the central frequency W£. 

After the discussion of the central spectral band, we now turn to the sideband frequencies around ll>l ± 2fi. Again, 
we consider detection in the far zone limit, but at this time the two detectors are at different positions. Under this 
conditions, the Cauchy-Schwarz parameters in Eq. (|I2p evaluate to 



XL = XR 



1 + cos(<5i - 8 2 ) 



-1 2 



3 + cos(<5i + 82) 



(14) 



This results is shown in Fig. 03a) for an interparticle distance r a {, = 5A and two atoms (N = 2). One may easily observe 
that for a large range of detector positions <5i,<52, the Cauchy-Schwarz inequalities are violated, i.e., xl = XR. < 1- 
These results for \l and xr can be easily understood by inspecting Eq. (fT2|) and Fig. [2] The Cauchy-Schwarz 



(2) ~* (2) ~~* 

parameters are given by the ratio of the product of the sideband second-order photon correlations g^KR) and gjift (R) 

to the cross correlation g LR {R) or g RL (R) squared. Therefore, both oscillatory structures in Fig. [5] are combined to 
give the result in Fig. 02a) . It is interesting to note that it is not possible to distinguish whether the first photon is 



0.1 0.2 0.3 0.1 0.2 0.3 

FIG. 4: (a) Second-order side-band photon correlations g^ c (R\, R2) and (b) g^ R (R\, R2) = jflfiJi,^) as function of the 
detector positions 01,02- The interatomic distance is r a b/\ = 5, and the number of atoms is N = 8. 

(2) — ' (2) ~* 

emitted on the left or on the right sideband, since the two different cross-correlations g LR (R) and g R [(R) are equal. 

V. MULTI-PARTICLE STRUCTURES 

In this final part, we discuss the case of N independent two-level emitters that are uniformly distributed in a regular 
chain. Thus the inter-particle separations between the neighbor particles is assumed uniform and equal to tq. As 
before, the atoms are driven by a strong laser field, ^> 1. We further assume detection in the far-zone limit, that 
is, the linear dimension of the chain L = (N — l)ro is much smaller than the distances between chain and detectors 
|i?i| and life |- Then the second-order correlation functions for the light emitted in the central- and the sidebands as 
well as to the cross-correlations can be evaluated to give 

^(i?i,ife) = l-jj + ^sWt+W + tfa-h)], (15a) 

si 2 2(ife,ife) = S$t(iMa) = l-^ + ^I^i-^), (15b) 

<?i 2 !(ife,ife) - 9 { rI{Rx,R2) = 1 + + 82) , (15c) 

J 2 ) (Z. J$„\ - J*) 



<?^(ife,ife) = g^' c {R 1 ,R 2 ) = 1 for X e {L, R} , (15d) 

where (f)(6) = sin 2 (N 6/2)/ sin 2 (5/2). 

In the multi-particle case, the Cauchy-Schwarz parameters Eq. (|12[) for the side-band photon correlations are given 
by the following expression 



XL = XR 



N 2 — 2N + cf>(5i -5 2 y 2 



N 2 + cf)(5i + 62) 



(16) 



Figure [3Jb) shows this Cauchy-Schwarz parameter for a linear chain of N = 8 atoms. In this figure, the Cauchy- 
Schwarz parameters assume values from below unity up to 4. In the following, we discuss two special cases to 
demonstrate this more clearly. First, consider a two-photon detector with 5i = 82 = 5. Then Eq. (fT6|) reduces to 



Xl = Xr= { n 2 + ^(26) ) (17) 

which for N ^> 1 and 5 = mr (n = 0, 1, • • • ) tends to Xl(r) —> 1 while for 5 = ir(l + 2n)/2 we have Xl(r) - > 4. 
If, for instance, 61(62) = and 6 2 (6i) = <5 with N ^> 1 then Xl(r) g° es to unity as Xl{r) — [1 — 1/A] 2 when 
5 = mr (n = 0, 1, • ■ • ) and to Xl(r) — [1 — 2/N] 2 for 5 = ir(l + 2n)/2. Thus, violation of CSI is likely to occur for 
moderate atomic structures. 
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Fig-BJa) shows the central band second order correlation function g^ c (R\, R2) for N — 8 atoms versus the detection 
angles a\,a2- Comparing this figure with Fig.[3]Jb), it is apparent that the dark regions in Fig.[3Jb) that correspond 
to a super-poissonian photon emission statistics for a pair of photons from the left and the right spectral sideband 
[see Fig.UJb)] are regions of high correlation in Fig. Ufa) as well. Thus, in these detection directions, super-poissonian 
statistics is observed in the central spectral band and the cross-correlations of the sideband photons. 

Finally we note that the properties of the correlation functions found here are quite different from those obtained in 
a Dicke-type sample [25|, 0, H3] ■ For example, Dicke-type samples do not exhibit the spatial dependence originating 
from the regular structure discussed here, since the Dicke model involves the assumption that all atoms interact with 
the electromagnetic fields with the same phase. In particular, the first-order correlation function in our present model 
is linear in the number of atoms N, whereas it is proportional to the number of atoms squared in the Dicke model. The 
unnormalized second-order correlation function is proportional to ./V 4 and does not depend on the detection direction 
in Dicke-type samples. For regular structures, the unnormalized second-order correlation function depends at most 
on TV 2 for particular detection positions. This detector position dependence arises from the geometrical phase factors 
exp(z/c ■ f ab ), and not from collectivity. Constructive interfere of these phase factors gives rise to the dependence on 
N 2 for particular observation directions. 

VI. SUMMARY AND CONCLUSION 

In summary, the intensity-intensity correlations of photons scattered by a strongly pumped linear atomic structure 
was investigated in detail. For this, the spectrum of the scattered light was separated into different spectral bands 
which we treated independently. In the central band at the driving laser field frequency, the resonant two-particle 
second-order interference fringes for two-photon detection have double resolution in the strong-field case as compared 
to the weak-field pattern. We found violation of the spatial Cauchy-Schwarz inequalities for photons emitted into the 
spectral sidebands, and that it is impossible to predict which photon is emitted first. Multiparticle atomic structures 
significantly enhance the intensity of the emitted correlated photons, and exhibit properties different from those in 
Dicke-type multiparticle samples. 
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